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Abstract 

We present a theory of valence holes as Luttinger spinor based qubits in p-doped self-assembled 
quantum dots within the 4-band k ■ p formalism. The two qubit levels are identified with the two 
chiralities of the doubly degenerate ground state. We show that single qubit operations can be 
implemented with static magnetic field applied along the z axis (growth direction) for a z operation 
and with magnetic field in the quantum dot plane, x direction, for a x operation. The coupling of 
two dots and hence the double qubit operations are shown to be sensitive to the orientation of the 
two quantum dots. For vertical qubit arrays, there exists an optimal qubit separation suitable for 
the voltage control of qubit-qubit interactions. 
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I. INTRODUCTION 



There is currently interest in using single electron spins in quantum dots (QDs) for quan- 
tum information processing (QIP) applications 1-6 . In III — V semiconductor quantum dots, 
such as GaAs, spins of conduction band electrons couple to nuclear spins of the host mate- 
rial and suffer decoherence due to hyperfine interactions. However, valence hole states are 
built from atomic p-type orbitals, which are expected to minimize the hyperfine interaction 
between the hole and surrounding nuclear spins 7 . The spin of a valence hole as a qubit is 
expected to have longer coherence time. A theory of heavy holes as qubits has been devel- 
oped in the single band approximation 7-9 . The justification for the use of the heavy hole 
single band approximation is that strain and QD confinement suppress the coupling between 
heavy hole (HH) and light hole (LH) bands. Under such an assumption, the qubit levels are 
defined by the J z = +3/2 and J z = —3/2 HH states. 

In such a simple model, the exchange coupling J needed to generate entanglement between 
spins of two holes is proportional to t 2 /U, where t is the tunneling strength of the hole and 
U is the on-site Coulomb energy. However, recent theoretical^ and experimental™ work 
show that in two vertically coupled disk-like QDs the spin orbit (SO) coupling between HH 
and LH bands changes the sign of the effective tunneling matrix element t as a function of 
inter-dot distance. If the inter-dot distance is smaller than a critical value, the hole state 
is mainly symmetric. However, if the inter-dot distance passes the critical value, the hole 
state is mainly antisymmetric. These results suggest that a description of an array of hole- 
based qubits requires taking into account the hole tunneling between dots with SO coupling 
between HH and LH bands. The theory of valence holes in quantum dots with strong SO 
coupling has been developed already^—. In this theory hole spin is strongly coupled to 
orbital motion and valence hole states are treated as Luttinger spinors. The two HH ground 
states are replaced by two Luttinger spinors-^ with different chiralities. 

In this work, we develop a theory of qubits based on the chirality of the valence hole-based 
Luttinger spinors. We are particularly interested in square-like self-assembled quantum 
dots grown on nanotemplates^ , which due to scalable architecture may lead to quantum 
information processing devices. We start with a single hole confined in an isolated QD. We 
explicitly define the qubit as two chirality states of Luttinger spinor and show how to perform 
single qubit operations. Our theory reproduces results obtained earlier by Kyrychenko and 
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Kossuth. Next, we investigate the tunneling of a valence hole in two orientations of two 
coupled square QDs as illustrated in Fig [TJ The tunneling barriers for vertically coupled 
quantum dots (VCQDs) and laterally coupled quantum dots (LCQDs) are modelled as finite 
potential wells, as shown in FigdJ The tunnel barrier strength is characterized by the band- 
offset between the dot and the barrier material. For VCQDs, we verify the results^ obtained 
for disk quantum dots showing the reversal of ground state character from symmetric to 
antisymmetric as a function of inter-dot distance. Close to the reversal, the vanishing 
of tunneling of a Luttinger spinor hole in VCQDs allows the benefit of strong confinement 
defined by growth with the possibility to control tunneling and hence double qubit operations 
using additional metallic gates. 



II. MODEL 



Following Ref. 
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we expand the wavefunction of a valence hole confined to a nanos- 
tructure defined by confining potential V(r) in terms of J z = 3/2,-1/2,1/2,-3/2 Bloch 
basis functions. The 4-band Luttinger Kohn (LK) Hamiltonian 19 reads, 

/ P+ R —S \ 

R* P- S 
-S* P- R 
\ S* R* P+ J 

where I is the identity matrix, fi is Planck's constant, Q = TieBjm^c is the cyclotron energy, 
and k is a material parameter. V(r) represents a 3D infinite potential well in case of a single 
QD in this study. The operator J is the angular momentum operator for spin 3/2 particle 
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where 71,2,3 represent the Luttinger material parameters, 723 = (72 + 7 3)/2, U a = k a — 
(e/ch) (A) a , k a = —i-§^-, a = x,y,z, and k- = k x — ik y . For this study, we choose In- 



GaAs/GaAs QDs. The InGaAs Luttinger material parameters are 71 = 11.01, 72 = 4.18, 
and 73 = 4.84. The subscript means no external field in this study. 

The LK Hamiltonian in Eq.flTJ exhibits several symmetries that we use to characterize 
the Luttinger spinors (the eigenstates of the LK Hamiltonian). The confining potential 
along the z direction is symmetric with respect to the origin defined in the middle of the 
structure; yet, due to spin-orbit coupling between different bands, our system does not have 
definite parity symmetry along the z direction but the time reversal symmetry of the system 
demands the two HH bands have opposite parity. The same also holds for LH bands. This 
allows us to define the chirality operator—^ which reads, 

/ i z \ 

2, 

Xz = . , (3) 

-i t 

\o -ij 

where i z is the inversion operator with respect to the z variable. The \ z operator has 2 
eigenvalues, which we denote by f, and \.. In the case of Xz =1\ the Luttinger spinors have 
even parity with respect to z for the first two components of the spinor and odd parity with 
respect to z for the last two components of the spinor. In the case of Xz =1, the parity pattern 
is reversed between the first two and last two components of the spinor. Furthermore, square 
QDs also have parity symmetry in the x—y plane, and we may define another parity operator 
Xxy analogous to the Xz by replacing i z with i x i y on the diagonal terms in Eq.Q where % x 
and i y are inversion operators along the x and the y direction respectively. Arguments for 
the Xz operator also apply to the Xxy operator. The only difference is that the R operator 
in Eq.(J2]) involves {kJ) 2 term which prevents further separation of x and y variables. Again, 
Xxy has 2 eigenvalues denoted by Xxy = +1 and Xxy = — 1- For Xxy = +1, the first two 
components of the spinor must have the same parity with respect to x and y, and the last two 
components of the spinor must have opposite parity with respect to x and y. For Xxy = — 1? 
the Luttinger spinors again have the parity pattern switched between the first two and last 
two components. 

The 4-band valence hole spinors in a square QD have well defined structures. The Bloch 
function part of the Luttinger spinors is determined by the lattice symmetries; while the 
envelope function part of Luttinger spinors is obtained from Eq.flTJ. In our analysis, we 
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expand the envelope function of the Luttinger spinors in eigenstates of the 3D infinite po- 
tential well of the size of our computational box. As mentioned, the chirality and the x — y 
parity symmetries allow us to label Luttinger spinors by Xz, Xxy, and N, which index the 
N th eigenstate of subspace Xz and Xxy F° r instance, the N th state of the subspace Xz =t 
and Xxy = +1 reads 

/ 
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where p = 1,2,3,... are positive integers. C*^ l f xy,N ' Jz represents the coefficient for the 
envelope function of a specific component \ J Z ) of Luttinger spinors. In our study, we assume 
that the confinement along the z direction is strong, so we only need to consider the first 
2 eigenfunctions of the ID infinite potential well in this direction. £ Sr (r) = */ cos (^p-^j 
if s r is odd and sin (^r^j if s r is even. W r is the computation box length along the 
r direction. \J Z = ±|,±|) represents the Bloch functions^ of zincblende structure with 
the total angular momentum projection J z . The Luttinger spinors in other subspaces have 
distinguishable parity patterns derived from Eq.Q by switching either the parity pattern 
with respect to the x — y variable or the z variable between first two and last two components 
of the spinor. 

From here onwards, the Luttinger spinor is written with the components corresponding 
to the Bloch functions in the order J z = 3/2, —1/2, +1/2, —3/2 as done in Eq.(j4j), and we 
will no longer explicitly write out the Bloch functions \ J Z ) as we have done above. 



III. HOLE STATES IN SINGLE QUANTUM DOT 

In this analysis, we take the entire Luttinger spinor as a representation of valence 
hole state confined in the quantum dots as opposed to adopting the HH single band 
approximation^. Due to the time reversal symmetry, we have a doubly degenerate ground 
state with distinct chirality at zero field. We encode the qubit with the chirality, Xz, of the 
ground state Luttinger spinors. We show that static magnetic fields applied along the z and 
the x direction indeed act as a x and a z operators respectively in the qubit subspace. This 
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result assures all the basic ingredients needed to perform arbitrary single qubit rotations. 



A. Hole under B z field 

By diagonalizing the LK Hamiltonian for a square QD, we find that the 2 ground states are 
characterized by quantum numbers (x z =t,Xxy = — 1) and (xz —\-iXxy = +1) respectively 
Let us denote the state with \\ z =T) by an d the other state with \x z =4) by |0). We 
give an approximate form of the envelope wave functions 



where = y cos \J%rJ and £2(7") = J -^r sin ^fp-J- |0) is obtained by applying the 
time reversal operator to |1). According to Kramer's degeneracy theorem, |0) will have the 
same set of envelope functions in Eq.(j5]) but complex conjugated and in a reversed order in 
the spinor. Due to the strong confinement along the z direction, the ground states of the 
system will contain a dominant contribution from the HH states. Therefore, the states |0) 
and |1) correspond to either a +3/2 HH or —3/2 HH occupying the (1, 1, 1) orbital of a 3D 
infinite potential box. 

To analyze the system in the presence of external fields B, we use the gauge A = 
{—B z y /2, B z x/2, B x y — B y x). First, we consider the system under an external field B z . 
Fig [2^ shows the energy spectrum of a single square dot charged with one valence hole as a 
function of B z . We see that the qubit subspace, the 2 lowest energy levels, are well isolated 
from the rest of the energy spectrum across a wide range of applied field. Due to the weak 
coupling between the qubit subspace and other excited states, we may treat the effects of B z 
fields by the first order Lowdin perturbation theory^, which is simply a projection of the full 
LK Hamiltonian into the qubit subspace. The entire LK Hamiltonian is decomposed into 
2 parts: unperturbed Hamiltonian at the zero field and additional perturbing Hamiltonian 
due to the field. The additional Hamiltonian with external field B z leads to ^ operators 
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of the form 

h 2 

Pi = ~ — (7i + 72) [(uj z ) 2 (x 2 + y 2 ) + u z (xk y - yk x )] , 

ZTTIq 

h 2 

Pi = 7, — (7i - 72) [(uzf(x 2 + y 2 ) + uJz(xk y - yk x )] , 

ZTTIq 

(6) 

where to z = ^JK Whether the applied field can break the chirality symmetry \ z and the x-y 
parity symmetry \ xy depends on the commutation relations between P^~^ operators and 
inversion operators i z and i x i y . Since P± commute with i z , the Hamiltonian can not mix 
the 2 spinors of different chiralities. Indeed, when we project the Hamiltonian we find the 
off-diagonal matrix elements (1\Hlk\Q) — 0, and the diagonal matrix elements (1\Hlk\^-) 
and (0\Hlk\Q) differ only by the Zeeman energy. The diagonal matrix elements read, 

(1|£l*|1> = ^(|a| 2 -|rf!| 2 -|d 2 | 2 )(3/2) + 
/d](-|6| 2 + | Cl | 2 + |c 2 | 2 )(l/2), 
(Q\H LK \0) = -(1\H LK \1), (7) 

where the coefficients a,b,c,d are defined in Eq.(jnj). As the ground state of flat QDs has 
dominant contributions from HH components, typically \a\ 2 0.9 is much greater than the 
magnitude of other coefficients. In HH single band approximation, one simply sets \a\ — 1 
and other coefficients to be 0. 



B. Hole under B x field 

Next, we consider the system subject to a constant B x field. Fig [2b shows the energy 
spectrum of a single square QD charged with one valence hole as a function of B x . The 
two qubit states are again well isolated from the rest of the energy spectrum across a wide 
range of B x . Hence, we repeat the same procedure to derive an effective Hamiltonian for the 
qubit. The vector potential is A = (0, 0, B x y) and the ' operators for the additional 
LK Hamiltonian with B x field read 

h 2 

Pi = 7^ — (71 - 272) [M 2 y 2 - 2u x yk z ] , 

2m 
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Pi = 7; — (71 + 272) [(wx) 2 y 2 - iu x yk z } , 
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where u x = We find that the diagonal matrix elements (1\Hlk\^) = (0\Hlk\Q), and 
the off-diagonal matrix elements are given by 

(i\h lk \o) = (^ 3/2 iA + l^ 3/2 ) + (^ /2 lA + l^ 3/2 ) + 
(^iPrb^ + ^iA-r 172 )- 

(g- l ' 2 W) + (g*' 2 \S*\g-^) - 
(^ 3/2 |5|^ 2 ) + (^ 2 |S|^ 3 / 2 ), 
w n x ( ll -2 l2 )^Im(d* 1 a). (9) 

where fl x = heB x /mc. \g l ) stands for the envelope function of the i hole states in Eq.Q, for 
instance, (r\g 3 ^ 2 ) = a^ 1 (x)^i(y)$ i i(z). The matrix elements (g~ 3 ^ 2 \P^~\g 3 ^ 2 ) do not suggest 
coupling of +3/2 HH and -3/2 HH. Rather, it actually represents the coupling of a chirality 
U P (Xz =T) +3/2 HH with a chirality down (x z =4) +3/2 HH. Due to time reversal sym- 
metry, the chirality down +3/2 HH must have the same envelope function as the chirality 
up -3/2 HH. Similar argument applies to all other matrix elements in Eq.Q. We further 
remark that the off-diagonal matrix element contains pairs of complex conjugates; thus, the 
matrix element is real-valued. 

The mixing of the qubit states is due to the simultaneous breaking of the parity symmetry 
in the x — y plane and the inversion symmetry along the z direction; for instance, we refer 
to the term yk z in both P x + and P{ operators in Eq.fjHJ) as responsible for the breaking of 
symmetries. The fact that we may couple the qubit states by static fields is in contradiction 
to analysis done in the HH single band approximation 8 . In HH single band approximation, 
the single qubit operation cannot be done by electron spin resonance (ESR) techniques be- 
cause the magnetic field cannot couple two HH states (in the leading dipole approximation). 
More sophisticated techniques such as electric dipole spin resonance (EDSR) are needed. 
However, in our proposal, we define qubit states with the entire Luttinger spinor, which 
is an admixture of all HHs and LHs, and the mixing of the two qubit states is due to the 
coupling between the two Luttinger spinors as manifested in Eq.Q. 

The a z operator for our proposed Luttinger spinor based qubit is essentially driven by the 
Zeeman energy; whereas the a x operator is based on the mixing of valence hole components 
by vector potential A and SO coupling. Due to the different mechanisms of how the 2 qubit 
states are operated on by B z and B x field, the effective g factor (for a simple spin in an 
external field) seems to have much stronger transverse component- 21 , g±_ » gn. Hence it 
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takes more time to perform a spin flip process with the Luttinger spinor based qubit. We will 
characterize the single qubit operating time with the h/AE gap , where AE gap is the energy- 
difference between ground and first excited state under B x field. From Fig [2b, we see the 
characteristic time is around 1.5 ns when B x = IT. However, modulating the strain of the 
host material can relax the QD confinement along the z direction, induce a stronger planar 
component of the g factors and improve the single qubit operating time. 

Finally, we remark that the Luttinger spinor based qubits are susceptible to the same 
:oherence that are already discussed for valence hole confined in QDs in 
However, electric field fluctuation in the background of the host material 
will not affect the Luttinger spinor based qubits. As the electric field induced dipole tran- 
sitions, (ipj\eE(t) ■ r\ipi), conserve the time reversal symmetry, the two qubit states, which 
are Kramers doublets, will not be mixed. 
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IV. HOLE STATES IN COUPLED QUANTUM DOTS 

We now turn to the analysis of double qubit operations with the Luttinger spinors based 
qubits. Our proposal for double qubit operations with the valence holes-based qubits relies on 
the following assumptions. First, we consider 2 valence holes localized in two different QDs 
in the weak tunneling limit and only the on-site Coulomb interaction is taken into account. 
Second, we consider that each QD only contains the two relevant qubit states. Under these 
assumptions, the hole-hole interacting Hamiltonian in second quantization reads^ 



H 2 h - ^2 e iP c )p c jp + t (°lp' c iP + c )p c i 

3,P 3 

+ 2 ^ipJ2pi3pi4p c jip C ] 2 p C i3p c i4p! (10) 
hv.hv 

J3pJ4p 

p 

where p are indices for QD number 1 and 2, p ^ p', ej p is the energy of valence hole state 
\j) on the p — th dot, t is the tunneling parameter between QDs, and Uj 1 j 2 j 3 j 4 is the on-site 
Coulomb interaction, with the matrix elements given in the appendix. This Hamiltonian 
can be greatly simplified because the Coulomb interaction conserves the chirality of 2 holes 
in the following manner: (a) if Xzi = Xz 2 and Xz z = Xz 4 , U jlhjsj4 = U c , (b) if Xz x ^ Xz 2 
and Xz 3 7^ Xzi, Uj 1 j 2 j 3 j i = U x , (c) if just one of the valence holes switches chirality, U = 0. 
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Taking also into account that we only consider 2 states of distinct chirality on each QD, the 
on-site Coulomb interactions term in Eq. ffTO]) reduces to^ 

U = -(Uc + U^hjfif, (11) 

where is the number operator, and j,j' represents the 2 chirality states on a QD. Differ- 
ent from electron spins, the on-site Coulomb interaction between valence holes is composed 
of a direct Coulomb term U c and an exchange term U x which enhance the overall interaction. 
However, apart from this difference, once we have parameterized the tunneling parameter t 
and U, we have the Hubbard Hamiltonian. In the strong Coulomb regime, the interaction 
between 2 localized particles in the Hubbard model can be reduced to two spins with the 
Heisenberg exchange constants J proportional to t 2 /U. 

An interesting phenomenon^ associated with valence holes is the possibility of engineering 
t to be either positive, zero or negative in a stack of vertically coupled cylindrical QDs by 
simply tuning the inter-dot distance. As the exchange coupling between 2 qubits is directly 
proportional to t 2 , this implies qubit-qubit interactions can be turned on and off as needed. 
We would like to investigate whether the similar phenomenon will happen with square QDs 
stacked either in a vertical or lateral structure. 

To estimate the magnitude of t, we rely on a tight-binding picture in which a valence 
hole tunnels between the two dots and the hybridization of local orbitals gives a symmetric 
state, l/\/2(|l) + |2)) with energy E — t and an antisymmetric state, 1/\^2(\1) — |2)) with 
energy E + t. Thus, the energy gap between first excited state and ground state is 2t. So we 
compute the energy spectrum of a single valence hole in a double QDs by exact diagonalizing 
Eq.((T]), then we extract the value of t from the calculated energy gap. 

A. Hole states in vertically coupled quantum dots 

The potential V(r) of VCQDs is modeled with an infinite potential well along the x and 
the y direction and double well potential profile along the z direction. In the barrier region 
between the dots, we set a constant finite potential of 320 meV corresponding to the band 
offset between InGaAs and GaAs, same as reported in Ref. |l(J, so we may draw a comparison 
between square-like and disk-like dots. 

As shown in Fig the lowest energy levels of different chirality subspaces cross at a 
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certain inter-dot distance. The insets of Fig [3^l present the most dominant HH wavefunction 
profile along the z direction of the ground state before and after the crossing point. The 
wavefunction plots show the reversal of symmetry for the ground state of the system. Similar 
phenomenon occurs in both square and disk dots. These results confirm that the effects of 
HH and LH mixing is insensitive to the confining potential in the x — y plane. These results 
also imply that we can shut down any undesirable hole - hole interaction between holes 
localized in different dots if the inter-dot distance is chosen at the point where the crossing 
occurs in Fig [3k. 

B. Hole states in laterally coupled quantum dots 

The potential of LCQDs is modeled similar to that of VCQDs except that the double 
well potential profile lies along the y direction. Our result, Fig [3b, shows no crossing of 
lowest energy states from distinct chirality subspaces. We understand that the confinement 
strength along the z direction plays a crucial role in determining how much HH and LH 
bands mixing a system will experience, because QDs are quasi-two dimensional devices with 
a much smaller dimension along the z direction. In the case of VCQDs, the coupled dots 
relax significantly the confinement strength of the hole states along the z direction and bring 
the LH and HH energies closer. However, a LCQDs structure relaxes confinement strength 
in the x — y plane which does not facilitate the mixing of HH and LH states. Hence no 
reversal of the ground state as a function of the distance is observed. 

V. CONCLUSION 

In summary, we consider the Luttinger spinor description of confined valence hole states 
in QDs. We identify the two qubit levels with two chiralities of the lowest energy Kramers 
doublet and suggest that B z and B x fields act analogously to the \ z and \x operators. 
For arrays of hole qubits we study tunneling of holes in vertical and lateral architecture as 
dominant mechanism for the qubit exchange interaction. We show that tunneling can be 
arrested for vertical pairs of quantum dots but not for lateral architecture. The capability 
to switch the sign of the effective tunneling t in VCQDs is demonstrated. This suggests the 
possibility of turning off the exchange interaction. With such exchange interaction being 
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very small one can envisage tuning this interaction with additional metallic gates. 
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Appendix: On-site Coulomb matrix elements 

The on-site Coulomb matrix elements among valence hole states \ji) = \XziXx yi , Ni) can 
be evaluated as follows 

2 

tj _ e V"^ s~l*Xz 1 ,Xxy 1 ,Nl,J zl „*Xz 2 ,Xxy 2 ,N 2 ,J Z2 „Xz 3 ,Xxy3,N3, Jz 2( ^Xz 4 , Xxy 4 ,N4,Jz 1 

31323334 ~ g n 3^ / j / j ni,mi,h n 2 ,m 2 ,l 2 ^113,1113,13 ^ 714,014^4 

Jz-\,Jzo 

1=1,2,3,4 

J^J^J^ ^ ^ g 2 "^ ^ nin2n 3™ 4 i^i^i Qx)G mirn2m3rri4 (y, y , Qy)Gi l i 2 i 3 i 4 (z, z , Qz), (A.l) 
where e is the static dielectric constants of InGaAs. We define 

G SlS2S3S4 (r,r',q r ) = jj drdr^U^U^iry^' 1 . (A.2) 
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FIG. 1: Schematics of (a) square vertically coupled quantum dots (VCQDs) of side lengths L=20nm 
and height H=2nm separated by a potential barrier Vq of width D v along the vertical direction, 
(b) Square laterally coupled quantum dots (LCQDs) of same side lengths L, height H separated 
by a potential barrier Vq of width along the horizontal direction. 
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FIG. 2: Energy levels of the 6 lowest lying hole states as a function of (a) magnetic field B z (applied 
in the z direction) and (b) magnetic field B x (applied in the x direction). 
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FIG. 3: Energy levels of the two lowest hole states as a function of the inter-dot distance in: (a) 
VCQDs (vertically coupled quantum dots) and (b) LCQDs (laterally coupled quantum dots). The 
critical distance at which the tunneling element is zero occurs around 1.8 nm in (a). The solid 
curve represents Xz =t state and the dashed curve represents the Xz =i state. Inset: The dominant 
heavy hole component of the ground state before and after the crossing. 
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